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A selection of recent results on heavy quark production at the HERA collider by the HI and ZEUS collaborations 
are presented with a focus on charmonium production in DIS, charm fragmentation and beauty production. 



1. Introduction 

In the heavy quark production processes there 
is at least one hard scale and hence perturbative 
QCD calculations are expected to be reliable. At 
HERA, where electrons or positrons of 27.5 GeV 
are coUided with protons of 920 GeV (820 GeV be- 
fore year 1998), pQCD predictions are tested in 
ep interactions. The HERA collider is a unique 
facility to investigate the parton dynamics, the 
heavy quark production and fragmentation prop- 
erties and other topics. 

The main kinematic variables are the virtuality 
of the exchanged photon, Q^, the photon-proton 
centre-of-mass energy, VF, the Bjorken scaling 
variable, x, and another scaling variable, y, which 
is the energy fraction of the lepton beam trans- 
ferred to the exchanged photon. 

Both the photoproduction {Q^ < IGeV^) and 
deep inelastic scattering (DIS, > IGeV^) 
regimes are studied. 

2. Inelastic J/ip production 

In Fig.^recent results on the inelastic J/ijj pro- 
duction in DIS regime measured by ZEUS T] are 
compared to different theoretical predictions in- 
cluding colour singlet (CS) and colour octet (CO) 
contributions as well as the CS contributions only. 
Neither shapes nor normalisation are described if 
CO contributions are included whereas CS contri- 
butions alone generally agree with the data. This 
is also true for the properties of the hadronic fi- 
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Figure 1. Absolute and normalised differential 
cross sections for the reaction ep — > eJ /^X in the 
kinematic region 2 < < 80GeV^, 50 < < 
250 GeV, 0.2 < z < 0.9 and -1.6 < Y^h < 1.3 as 
a function of log(M|-) and Yx- The data are com- 
pared to LO NRQCD (CS-fCO) predictions, a 
LO NRQCD (CS) calculation, a prediction in the 
fc^-factorisation approach within the GSM and 
the Cascade (CCFM-based) MC predictions. 

nal state (invariant mass, Afx, and rapidity, Yx) 
which were measured in this process for the first 
time. From similar studies in photoproduction |2] 
it is known that rcsummation techniques for soft 
gluon emission can improve the theoretical de- 
scription, however such calculations are not avail- 
able in DIS. Cascade Monte Carlo predictions 
are above the data but the shapes are well de- 
scribed. A calculation in the fcx-factorisation ap- 
proach based on BFKL evolution equations gives 
the best description of the data. 



2 



I. Katkov 



3. Charm fragmentation 

Fragmentation ratios and fractions in photo- 
production has been analysed by HI 3 and 
ZEUS 4^. In this analysis the charm gromid 
states, D+, A+, and the charm vec- 

tor meson D*+ were measured. As an ex- 
ample in Fig. 121 the three-track invariant mass 
M{KpTr) distribution is shown for the A+/A~ 
candidates after subtraction of reflections from 
£)+ and decays. In charm fragmentation, 
the ratio of neutral to charged D-meson rates, 
Ru/d = {cu)/{cd) = 1.100±0.078(stat.), is consis- 
tent with the isospin invariance, the strangeness- 
suppression factor is 7^ = (2cs)/{cd -f cu) = 
0.257 ± 0.024(stat.) and the fraction of charged 
D mesons produced in a vector state is Py = 
0.566 ± 0.025(stat.) which is not consistent with 
the naive spin counting 4 . The fragmentation 
fractions, /(c — > I?, Ac), were also measured 0]: 
0.217 ± 0.014 for D+, 0.523 ± 0.021 for D", 
0.095 ± 0.008 for D+ , 0.144 ± 0.022 for A+ and 
0.200 ± 0.009 for D*+ 2. These values are consis- 
tent with HI results '3'. Also they are in general 
agreement and have precision competitive with 
the combined e~^e~ data ;5' thus supporting the 
universality of c-fragmentation. 

A study of the charm fragmentation function 
have been done by HI in DIS (2 < < 
100 GeV^ 0.05 < y < 0.7). In the definition 
of the fractional transfer of the c quark energy 
to a D* meson {D* Kirirs, pT > 1.5 GeV, 
jr/l < 1.5) the c quark energy was approximated 
by the energy in the appropriately defined hemi- 
sphere associated to the D* in the j*p frame: 
Zhem = [E + Pl)d' /Y^hemi^ + Pl)- This was 
done to establish a close analogy to the mea- 
surements in e+e~ experiments. The distribution 
over z is in agreement with recent e+e~ annihila- 
tion data from CLEO at a similar centre-of-mass 
energy of the cc pair. The data are well described 
by Rapgap Monte Carlo where either Peterson 
or Kartvelishvili parametrisation of the fragmen- 
tation function was implemented. The best de- 
scription was achieved with e = 0.018to^oo4 (Pe- 
terson) or a = 5.9lQ'g (Kartvelishvili). The e 



value is lower than the previous ZEUS photopro- 
duction result e = 0.064 ± 0.006l°;°J^ 0. 
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ratios and fractions; all others are omitted. 



Figure 2. The M{KpTr) distribution for the 
A+ / A~ candidates in photoproduction. 



4. Charm with jets 

Heavy quark production in association with 
jets measured recently by HI and ZEUS [2] 
provides a test of a multiscale process. The ZEUS 
measurement of jet cross sections in charm photo- 
production (Q2 < lGeV^ 130 < ly < 280 GeV) 
has been confronted with both massive and mass- 
less calculations. In the so-called massive scheme 
there are three active quark flavours in the pro- 
ton or photon with heavy quarks produced dy- 
namically via e.g. boson gluon fusion. The ap- 
proach is expected to be valid at low transverse 
momenta of the c quark. In the resummed mass- 
less scheme there are up to five active partons and 
it is valid at high transverse momenta. Combined 
calculation schemes (like variable flavour number 
scheme) also exist. 

The D* mesons {D* —>■ Kirns) were required to 
be in the kinematic region px > 3 GeV, \ri\ < 1.5. 
Jets (fcy cluster algorithm) with —1.5 < 7yj°* < 
2.4, £;j^* > 6 GeV were selected. For the dijet 
analysis the highest Ei^*" jet was required in ad- 
dition to satisfy Ei^^ > 7 GeV. 

For all Ei^^ and tf'^^ the shapes of inclusive jet 
cross sections with and without a D* tag are well 
described by both theoretical predictions whereas 
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Figure 3. Dijet correlations in charm photo- 
production compared to massive NLO QCD pre- 
dictions. The beauty component in the data 
(Herwig MC estimate) is also shown. 

the normalisation is described by their upper lim- 
its. There is no excess in the forward region (pos- 
itive rf°^) observed previously in the purely in- 
clusive jet measurements. In Fig. |21 dijet corre- 
lations are shown in comparison to the massive 
NLO QCD predictions: the fraction of the pho- 
ton momentum participating in the dijet produc- 
tion, x°^^; the azimuthal separation of the jets, 
A(j!)JJ; the transverse momentum squared of the 
jets, (p^^)^; and the dijet invariant mass, A'PK 
The dijet correlations are directly sensitive to the 
higher order corrections as for the lowest order 
2^2 process Acj^^ = tt and {p'j')'^ — 0. The mas- 
sive predictions are low in comparison to the data 
at low A(/)jj and high (pJ^)^. This discrepancy is 
enhanced for the resolved-enriched {x'^^ < 0.75) 
sample. The Herwig MC describes the shape of 
the measurements well. This indicates that ei- 
ther even higher-order calculations or a matching 
of the NLO matrix elements to a parton shower 
program (MC@NLO) are needed. 

Jet shape variables in charm photoproduction 
(g2 < lGeV^ 0.2 < y < 0.8) have been anal- 
ysed by HI ilOj. An insight into the hard scat- 
ter process is possible due to the fact that the 
integrated jet shape ^p{r) is expected to be differ- 
ent for quark- and gluon-initiated jets. Jets were 
reconstructed with the kx cluster algorithm. In 
the dijet sample (p^'"''"' > 7(6) GeV, < 1.7) 
one jet was tagged by a semileptonic decay muon 



(p^ > 2.5 GeV) as being initiated by a c quark 
using the two-dimensional fit of the muon trans- 
verse momentum with respect to the closest jet, 
p^'' , and the muon impact parameter, 6, dis- 
tributions to a mixture of ones in Pythia MC 
for light, charm and beauty quarks. The cut 
p^'' < 1 GeV was imposed to get a charm en- 
riched event sample. In Fig. 0] the average in- 
tegrated jet shapes {ip{r)) are shown for the jet 
without the charm tag. For the resolved enriched 
sample {x°^^ < 0.75) gluon jets are expected to 
dominate via the resolved photon charm excita- 
tion process, g ^ eg, according to the Pythia 
MC predictions. This dominance is not observed 
in the data and on the contrary the integrated 
jet shape {ip{r)) is similar to one in the direct 
enriched sample (a;°'^'^ > 0.75) where quark jets 
are expected to dominate via the photon gluon 
process, 75 ^ cc. 
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Figure 4. Average integrated jet shape for jets 
without charm tag for two different regions of 
^obs (resolved- and direct-enriched). The data 
are compared to Pythia MC predictions (total, 
direct and resolved contributions). 



5. Ff' and Ff 

The inclusive lifetime tagging technique have 
been used by HI for the identification of heavy 
quark hadrons. The method is based on the track 
impact parameter, (S, as measured by the HI ver- 
tex detector. For heavy quarks tracks have high 
significance, S, defined as (5 divided by its error. 
Hence significance distributions Si are considered 
where Si is the significance of track with i-th high- 
est absolute significance in an event. The c, h 
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and light quark fractions are extracted from the 
simultaneous fit to at least two distributions Si 
and 5*2. The method features low extrapolation 
factors (high acceptances). 



Q-= 3.5 GcV- 
• HI Preliminary 
▼ III Dala (High Q^) - 
--- MRST04 

CTIiQ6HQ - 
— CCKM - 


12 GcV^ 


-_ Q^= 25 GeV- J 
^ 1 1 1 1 iH^^^i-V^ 


■_ ..Q-= 60 GeV^ J 
^ I^^'^^^M^^ 


Q^= 200 GcV- ~ 

- X: - ' 


Q-= 65.0 GcV" " 



Figure 5. The structure function F2 as a func- 
tion of X for different values of Q^. The data are 
compared to NLO predictions in the framework of 
the VFNS (MRST, CTEQ) and predictions based 
on CCFM parton evolution. 

The method allowed a purely inclusive mea- 
surement of charm and beauty dijet cross sec- 
tions ^2 as well as the first measurement of 
the beauty contribution in the structure function. 
Using three distributions Si, S2 and 5*3, the mea- 
surement of FI'' shown in Fig. was extended to 
as low as ^ 60GeV^ [12 . The data are 
described by the variable flavour number scheme 
(VFNS) NLO calculations. In the kinematic re- 
gion of the measurement the beauty cross section 
is on average 0.8% of the total ep cross section 
(2.7% in the high-Q^ measurement). 

6. Beauty with dimuons 

Specific topologies of events with B decays al- 
low the identification of heavy quarks. Following 
the strategy of the previously done fi—D* cor- 



relation analyses, the dimuon charge correlations 
and invariant mass distributions have been stud- 
ied by ZEUS to measure beauty production jl3j . 
Despite rather restricted statistics the analysis 
is sensitive to B mesons almost at rest due to 
low background and high acceptance down to 
low muon transverse momentum, pt^ > 1.5 GeV. 
The muon rapidity range, —2.2 < 77^ < 2.5, 
is also wide due to almost full rapidity cover- 
age of the muon detectors used in the analy- 
sis. The azimuthal separation of the the muons, 
A0^^, was measured for high dimuon invariant 
mass, m'^'^ > 3.25 GeV, thus providing a direct 
probe of the bb correlations. The total cross sec- 
tion measurement, crtotisp — > bbX) — 16.1 ± 
1.8(stat.)t4 g(syst.) nb, is to be compared to a 
NLO QCD prediction, 6.8t?;?nb. 

7. Summary and outlook 

The large data sample accumulated during 
HERA I running allowed HI and ZEUS to test 
thoroughly different aspects of QCD in charm 
and beauty production processes. The NLO QCD 
calculations, where available, were found to be 
in general agreement with the data, although 
some measurements suggest that even more pre- 
cision is needed. New results from the post- 
upgrade HERA II running with increased colli- 
sion rates, polarised beams and upgraded detec- 
tors will bring more experimental precision in the 
heavy flavour studies. 
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